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We study the collisional energy loss mechanism for particles produced off mass-shell in a finite 
size QCD medium. The off mass-shell effects introduced are to consider particles produced in wave 
packets instead of plane waves and the length scale associated to an in-medium particles' life-time. 
We show that these effects reduce the energy loss as compared to the case when the particles are 
described as freely propagating from the source. The reduction of the energy loss is stronger as this 
scale becomes of the order or smaller than the medium size. We discuss possible consequences of 
the result on the description of the energy loss process in the parton recombination scenario. 

PACS numbers: 12.38.Mh, 24.85.+p, 25.75.-q 
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I. INTRODUCTION 

The problem of collisional parton energy loss in a QCD 
medium has been revived due to recent RHIC data on 
non-photonic single electrons [l[ that are not well de- 
scribed within radiative energy loss calculations. Colli- 
sional energy loss has been a subject of research from long 
a g° @> S 0] • During the pioneering years, an important 
question was to understand how to handle the infrared 
singularities in perturbative calculations at finite temper- 
ature. The advent of resummation techniques clarified 
this point and allowed to reliably compute the energy 
loss of a heavy parton traversing an infinie medium to 
lowest order in perturbation theory [5|, |6fl. Short after, it 
was estimated that radiative energy-loss in a finite size 
medium was a more important mechanism to account for 
energy losses of energetic partons @, H, Nonetheless, 
even more recent studies [l(| El suggested that for a 
range of parameters relevant to RHIC energies, radiative 
and collisional energy losses for heavy quarks are of the 
same order of magnitude. These last calculations where 
done for infinite QCD media. The outstanding question 
was whether collisional energy loss for finite size media 
was also significant. 

In this context there where two results seemingly in 
contradiction [12, EH- In Ref. [l2} a semiclassical ap- 
proach based on linear response theory computes the col- 
lisional energy loss by means of the work done by the 
response chromoelectric field on the color-charged heavy 
parton traversing the medium. The infinite medium limit 
of this description agrees with the collisional energy loss 
result at high temperature -up to color factors- obtained 
from a perturbative approach using HTL effective propa- 
gators jjj . The original claim that the finite size medium 
induced energy loss is strongly suppressed compared to 
the infinite medium case was later revised by properly 
subtracting the kinetic energy associated to producing 
the particle within the medium (l9| . 



On the other hand, in Ref. [131 ] a lowest order perturba- 
tive calculation using HTL propagators finds that finite 
size effects on the collisional energy loss arc not signif- 
icantly suppressed as compared to the infinite medium 
case. The formulation of the problem is based on the 
assumption that the scattered particle originates within 
the medium but otherwise is produced on mass-shell. 

However, when particles arc emitted by sources lasting 
a finite amount of time they are not necessarily produced 
on their mass-shell since the source emits over a (wide) 
range of energies. A physical consequence is the possibil- 
ity that the particle losses its identity within the medium. 
In the midsts of a high energy heavy-ion collision, such 
possibility can be realized in the recombination of a jet 
parton with the partons from the surrounding medium. 
Recall that during the propagation inside a deconfincd 
QCD plasma, a fast parton can have not only induced 
gluon radiation but also induced absorption from ther- 
mal gluons. This process can fairly well be considered 
as parton recombination, which is one of the accepted 
mechanisms that are used to describe the distinct fea- 
tures of meson and baryon spectra that include a baryon 
to meson ratio larger than one for p t > 2 GeV in cen- 
tral Au + Au collisions at RHIC [lj] and their different 
azimuthal anisotropies jlEj . This point is addressed in 
the context of the modification of parton fragmentation 
functions induced by medium effects in Ref. [161 ]. When 
a parton recombines it certainly losses its original iden- 
tity and the energy loss should stop being described in 
terms of parton degrees of freedom. Parton recombina- 
tion from a jet with thermal partons to form intermediate 
Pt hadrons is a viable scenario in the case of light flavors, 
given the features of the proton to pion ratio and even for 
s quarks, given the features of the A to kaon ratio [l7l |. 
Other off mass-shell effects can be of relevance as well 
when studying if and how the virtuality of the propagat- 
ing parton affects the in- medium splitting functions [1 81 ] . 

In this work we study one of such off mass-shell ef- 
fects, namely a possible finite life-time of the scattered 
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partons -originating within and traversing the finite size 
medium- in the description of the collisional energy loss 
mechanism. We introduce the possibility that the scat- 
tered parton be described in terms of a propagator con- 
taining a parameter associated to the parton's life time. 
We argue that for recombining partons, this picture could 
be used to consider collisional energy losses only up to 
times when these recombine with thermal partons from 
the medium from where the energy loss process should 
start being described in terms of hadronic degrees of free- 
dom. 

The work is organized as follows: In Sec.|Il]we rewrite 
the expression for the collisional energy loss in a finite size 
QCD medium based on the formalism used in Ref. [13 |. 
allowing for particles being emitted off mass-shell and 
having a finite life-time. In Sec. lIIII we give the numerical 
estimates for the collisional energy losses of heavy and 
light quarks using parameters relevant for RHIC energies. 
We use the cases studied in Ref. [l3[ as a base to compare 
our results to. We finally conclude and give an outlook 
of the consequences of the result in Sec. HVl 



II. ENERGY LOSS 
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FIG. 1: Scattering diagram describing the lowest order con- 
tribution to the energy loss due to the collision of fermions 
with incoming and outgoing momenta P and P' with medium 
fermions with initial and final momenta K and K' , respec- 
tively. The blob in the intermediate gluon line represents the 
effective HTL propagator. Incoming fermions are produced 
by the source in wave packets and not necessarily on their 
mass-shell. 



We start by describing the elemental interaction of a 
fermion with momentum P M = (po, p) (not necessarily on 
its mass-shell), magnitude of its velocity v = |v| = p/E, 
mass M and spin s with a massless fermion in the non- 
expanding medium, with momentum = (fc,k) and 
spin A, through the exchange of a gauge boson with mo- 
mentum = (w,q) and by means of a coupling con- 
stant g. For elastic collisions, these particles retain their 
identities and after the scattering they have momenta 
P ni = (E',p') and K" 1 = (fc',k') and spins s' and A', 
respectively, and E' = \J p' 2 + M 2 . The scattering dia- 
gram associated to the process is depicted in fig.[TJ When 
the incoming massive fermion is on its mass-shell, i.e. 
Po = E = ^/ p 2 + M 2 , the expression for the matrix ele- 
ment describing this process is given by (see Eq. (6) in 
Ref. 0) 



iM = ~g 2 J d 4 x j(t,x.)e lP x J d 4 Xl J d 4 x 2 
f d 3 p r d 4 q n ( a \ e "i-(xi-x2) 

x u(k', X')e iK '- X2 ^u{k, \)e- lK - X2 

x eih-t) 9(L/v-{h-t)), (1) 

where D a p will become the effective HTL gluon propaga- 
tor, j(t, x)e x is the amplitude associated to the source 
to produce an incoming particle with momentum P and 



x» = (i,x), 

x^ (t i , ) , i 1,2. 



The step functions in Eq. |1} represent the conditions 
that the produced fermion interacts within the plasma 
at ti after being produced at t and before it leaves the 
plasma at L/v + t, under the approximation that its ve- 
locity remains constant. 

When the source produces particles off their mass- 
shell, the matrix element describing the process can be 
written as 



iM = Jd 4 xj(t, x ) Jd 4 Xl Jd 

d 4 p f d 4 q 



•''2 

D a(n)p iq ' < ' Xl ~ X2 ^ 
(2tt) 4 J (2Tr)4 Uaf3(q)e 

x u(p', s')e lP '- Xl 1 a S(P)e lP - {x ^ Xl) u{p, s) 

x u(k' , X'yK'-^^uik, \)e~ tK - X2 

x 9(tx-t) 9(L/v - (ti - t)), (3) 

where S(P)e lP '( x ~ Xl " l u(p, s) represents the amplitude to 
propagate a fermion mode with momentum P from x to 
X\. Notice that Eq. © reduces to Eq. (TJ when S(P) -> 
S (P) = {2tt)(2M)K+{P)8{P 2 - M 2 )6( Po ), where 



A+(P) = 2«(p,*)«(p,s) 



M 



2M 



(4) 



(2) 



is the projector for positive-energy solutions. 

Let us however consider the situation where S(P) does 
not describe free, on-mass shell propagation, but instead 
the propagation of a wave packet with a finite width rj. 
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In order to explore a simple scenario, let us recall that 

^ ^c 1 _\£/o2 ~\/f2\ 



lim — , 

r,^0+ (po-E p ) 2 + if 



(2E p tt)S(P 2 - M 2 )9(p ), (5) 



where E p — +\/p 2 + M 2 . To consider a finite width, we 
take r\ finite and write 



S(P) = 



_ /2M 



V E p J ip - E p ) 2 + n 



rA+(P). 



(6) 



Upon the change of variable y = x\ — x and after inte- 
gration over d 4 X2, d 4 y, d 3 p, d 4 q and d 4 x in Eq. © we 
get 



2i[ — ) / dP0_ e -t(po-E'-u)L/2v 



' M '■'EpJJ_ 00 (2n) 

sm[(p - £" - cj)L/2w] 
X (p - E> - u;){(p ~ E p ) 2 + V 2 } 
x iM (p )j(P), (7) 

where j is the Fourier transform of j(t, x) and 

*-Mo(:Po) = ig 2 D aP {K'-K) 

x u(p',s')7 Q wb,s)w(fc',A')7 /3 w(fc,A) (8) 

and we used that A + (P)u(p, s) = u(p,s). -Mo(Po) is 
the matrix element describing the scattering process in 
infinite volume and for the following discussion, we have 
emphasized its dependence on p a . 

In order to perform the integral in Eq. j7]), we write 

1 1 



(po - E p ) 2 + V 2 [{p - E p ) + irj\ [(po - E p ) - ir)] ' 

and for convergence, close the contour of integration on 
the lower po complex half-plane, which picks the pole at 
Po = Ep—irj. This results in the expression for the matrix 
element 



iM 



e -r)L/2v e -i(E p —E' —u)L/2v 



sm[(E p —E'—u)- ir))L/2v] 
X (E p - E' - oj - iri) 

x iM (E p -ir))j{P). 



(10) 



Notice that for a consistent description, we require the 
condition rj <C E p meaning that the central energy of the 
wave packet is much larger than it's width. Therfore in 
Eq. (|10[) we can approximate 



M {E p -ir,)~M {E p ). 



(11) 



This approximation cannot be made for the rest of the 
factors in Eq. (flTJ)) since the term E p — E' — u> is of order 
of the transfered momentum that in turn is of order of 
the medium's temperature which may not be much larger 
than the wave packet's width. 



The square of the matrix element given in Eq. (|10[) . 
averaged over the initial spin s and summed over all other 
spins is therefore 



.A, A' 



-nL/v ( J_ 



sin[(w v ■ q + irf)L/2v] 



lo — v • q + ir\ 



\ Y, \Mo(E p )\ 2 , (12) 



,A,A' 



where we use that for an energetic incoming fermion 
E p — E' ~ v • q. It is worth mentioning that if instead of 
using the propagator in Eq. ©, one uses the free Feyn- 
man propagator, the result for the square of the matrix 
element, averaged over the initial and summed over all 
other spins, yields the result found in Ref. [l3j |. 

Hereafter, we specialize to the description of the scat- 
tering of the fermion (quark) in the QCD plasma. The 
differential energy loss dE is related to the differential 
collisional interaction rate dT by dE — uodT On the 
other hand, dT is given by 



d 3 NdT = ~ £ \M\ 2 



d 3 p' d 3 k d 3 k' 
(2n) 3 2E' (2n) 3 2k (27r) 3 2fc' 



s,s' ,A,A' 

x ni q {k)[l±ni q {k')] 

1 ,..,2 d 3 p' d 3 k d 3 k' 



2 s ^ ; '" " {2it) 3 2E> {2n) 3 2k (27r) 3 2fc' 



(13) 



where 

d 3 N 



->' . 2 , ' /, , .' > 2 



d R \j(P')\ 2 [u(p>, S ')u(p',s')Y 



{2ir) 3 2E> 



d R \j(P 



/ \ 1 2 



d 3 p' 



(2ir) 3 2E'' 



(14) 



represents the number of (non-scattered) particles into 
the phase space volume in the interval p' and p' + d 3 p' 
with d;?, = 3 for the SU(3) fundamental representation. 
Also, in Eq. (TT3]) we have used that when describing the 
collisional energy loss, n| ? (fc)[l ± nL(fc')] = n| (fe) since 
the term proportional to nf q (k)n^ q (k') is odd under the 
exchange of k and k' and integrates to zero [5( , and have 
defined 



E 



•eq 



(k). 



(15) 



S=q, q, g 



For a source producing energetic particles with a large 
spread in momentum, we can take the approximation 
|j(P')| 2 ~ |_7 (i- 3 ) | 2 . Therefore, the collisional energy loss, 
considering the finite width of the scattered wave packet 
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FIG. 2: (Color online) Fractional energy loss for (a) light, (b) 
charm and (c) bottom quarks as a function of their momenta 
for a fixed medium's length L — 5 fm. The uppermost curve 
in each case corresponds to the description without off mass- 
shell effects. This is compared to the case with off mass-shell 
effects for two values of r\ = 40, 80 MeV. In each case, the 
fractional energy loss decreases as the value of r\ increases. 



can be written as 

-r/L/v 

AE ~ C R - 



d 3 k 



El J (27r) 3 2fc 



(A:) 



d 3 k' 
(27r) 3 2fc'' 



sin[(w v q + irj)L/2v] 



U) — v • q + vq 
\ \Mo(E P )\ 2 . 



(16) 



Notice that Eq. ([16]) is modified with respect to the cor- 
responding expression in Ref. [13( by the ^-dependent 
exponential factor and the i] dependence in the argu- 
ments of the sine function and in the energy denomina- 
tor. When rj — > 0, the corresponding expression for the 
energy loss in Ref. [HI is recovered. In order to find the 
explicit expression for Eq. (fTT)]) . recall that the effective 
gluon propagator can be written as 



D» v = -P^At - QP^A L 



(17) 



where in the HTL approximation, the effective transverse 
and longitudinal gluon propagators are given by 



A" 1 



to 



2q 



JL 
2 

U3 - 



(uj 2 — q 2 )m 



2q 2 



uj + q 

LU 

+ 2i 



■In 



(18) 



where m 2 D — g 2 T 2 (l + Nf/6) is the square of the De- 
bye mass and, working in Coulomb gauge, the only non- 
vanishing components of the transverse and longitudinal 
projectors are 



pij — §ij 

Q 00 = l. 



q l q> 



(19) 



Using Eqs. (17]) -(]l9]).the averaged over initial, summed 
over all other spins, matrix element squared describing 
the underlying scattering process in vacuum is given by 

\ £ |M)| 2 = 16^j|A L ( 9 )| 2 (^' + k.k') 

s,s' .X,X' \ 

+ 2 Re [A(q) L A(q)?p] 

, v-qk'-q 



k v • k' 



k' \y • k 

|At( 9 )| 2 



r 

v • q k • q 
2 [ v • k • 



v • q k ■ q 



{kk' - k • k') ( v 



2 v q v q 



(20) 

Since for a non-expanding medium, the energy loss does 
not depend on the direction of v, we can simplify Eq. (flfl]) 
by averaging over the direction of v Q . This is most con- 
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(a) Ligh; 




the appendix. After averaging over the directions of v, 
the expression for the energy loss can be written as 



AE = <^El e -vL/v 

2n 4 jo 

( f k r 

x / qdq / udu) 



'eq 



{k)dk 



9max rq 1 

qdq / Loduo 

I. J q-lk , 



x |A L (g)| 



,{2k + uj) 2 -q 2 



Jo 



IAtI 2 



[q 2 - Lu 2 ][(2k + lu) 2 + q 2 } 



(b) Charm 




L (frr) 



x [{v 2 q 2 - lu 2 )J + 2uuJ x - J 2 ] 



(22) 



where the limits of integration over u> and q take into 
account that for the considered scattering amplitude, the 
transferred four-momentum is space-like and 



2k(l + k/E p ) 
l-v + 2k/E p ' 



(23) 



is obtained from the approximation that the maximum 
energy transfered occurs for backward scattering (2lj . 




L (fm) 



FIG. 3: (Color online) Fractional energy loss for (a) light, (b) 
charm and (c) bottom quarks as a function of the medium's 
length for a fixed quark momentum p = 10 GeV. The up- 
permost curve in each case corresponds to the description 
without off mass-shell effects. This is compared to the case 
with off mass-shell effects for two values of rj — 40, 80 MeV. 
In each case, the fractional energy loss decreases as the value 
of rj increases. 



veniently performed by introducing the auxiliary func- 
tions 



Ji 



dn 

47T 



sin[(w — v ■ q + irj)L/2v] 



lu — v • q + IT] 



(w-vq)*, (21) 



i = 1, 2, 3, in terms of which, the average over the 
different powers of v appearing in Eq. (|2T)|) can be ex- 
pressed. The functions in Eq. (|21[) are explicitly given in 



III. NUMERICAL RESULTS 

In order to present the quantitative behavior for the 
energy loss, we take standard values for the parameters 
involved. We give examples of the effect for light as well 
as for heavy flavors both, to study the mass effect and to 
directly compare to the findings of Ref . [l3| . The plasma 
temperature is taken as T = 0.225 GeV, the effective 
number of flavors Nt = 2.5, the strength of the coupling 
constant a — g 2 /47r = 0.3 and the Debye mass mr> = 
0.5 GeV. The bottom quark mass is taken as 4.5 GeV 
whereas the charm quark mass is taken as 1.2 GeV. We 
take the mass of the light quarks as 0.2 GeV. 

Figure [2] shows the fractional energy loss AE/E for 
light, charm and bottom quarks in a finite size medium 
with L = 5 fm, for the cases with and without off mass- 
shell effects. For the curves describing the off mass-shell 
effects, we consider two values, rj = 40 MeV and 77 = 80 
MeV. Notice that in all cases, a finite value of 77 produces 
the energy loss to decrease as compared to the case rj = 0. 
The decrease is more important for larger values of rj. 

Figure [3] shows the fractional energy loss for light, 
charm and bottom quarks as a function of the medium's 
length L for a fixed quark momentum p = 10 GeV, com- 
paring also the cases with and without off mass-shell ef- 
fects. For the curves describing the off mass-shell effects, 
we consider once more the two values, 77 = 40 MeV and 
7/ = 80 MeV. Notice that a finite value of rj causes the 
fractional energy loss to asymptotically reach a maxi- 
mum value as the medium's size increases. This is in 
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sharp contrast with the case where no off mass-shelf ef- 
fects are considered where for large L the energy loss 
increases linearly. 

Figure 0] shows the behavior of the fractional energy 
loss as a function of r\ also for light, charm and bottom 
quarks for a fixed value of the quark momentum p = 10 
GeV and a fixed value of the medium's size L = 5 fm. 
The fractional energy loss decreases with increasing rj and 
the decrease is similar in shape, regardless of the quark 
mass. 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have studied the off mass-shell effects 
on the collisional energy loss of particles, produced and 
scattered within a finite size QCD medium, associated 
to the introduction of a finite-width wave packet and 
therefore a finite particle's life-time. We have shown 
that this effect decreases the energy loss as compared to 
the case when these particles are produced on mass-shell 
and therefore live longer than the medium, fragmenting 
outside it. We have argued that this picture should be 
applied in particular to energetic partons that recom- 
bine with thermal partons and thus hadronize within the 
medium. 

Recall that the length scales playing a role for the 
energy loss mechanisms in a finite size, thermal, non- 
expanding medium are the medium's size L, the average 
distance between collisions d ~ l/T, the Debye radius 
rr>, the mean free path 5 ~ l/g 2 T and the particle's 
formation time tf ~ 1/E p . When considering the in- 
medium particle's lifetime, the length scale r\~ x is intro- 
duced. For the description of the scattering process in 
terms of a perturbative picture, it is required that the 
hierarchy of scales 

l/T < r D < S (24) 

be satisfied. The requirement that the scattering particle 
can be described in terms of an oscillating mode means 
that this mode is not damped too fast, which in turn 
translates into the condition 

1/E p « 1/r,, (25) 

and can be thought of as the condition that the time the 
medium takes to produce the particle is much shorter 
than the in-medium particle's lifetime. For very ener- 
getic partons, it is safe to assume that tf ~ 1/E p is the 
smallest of all length scales. Moreover, for media sizes of 
the order expected to be produced in relativistic heavy- 
ion collisions, it is also safe to assume that L is larger 
than the mean free path and therefore the hierarchy of 
scales 

1/E p < l/T < r D < S < L, (26) 

follows. In fact, the results in Refs. can be viewed 

as meaning that as long as the medium size is larger than 
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FIG. 4: (Color online) Fractional collisional energy loss for 
light, charm and bottom quarks as a function of r\ for a fixed 
quark momentum p = 10 GeV and a fixed medium's size 
L — 5 fm. The fractional energy loss decreases with increasing 
■q and the decrease is similar in shape independent of the quark 
mass. 



the Debye radius, the collisional energy loss for on mass- 
shell particles in a finite size medium is not suppressed 
as compared to the infinite medium case. 

The situation changes when introducing the in-medium 
particles' life-time. The results of this work show that 
when f? -1 < L the effects are strong. They cease to 
matter for r/ — > 0. 

It should be emphasized that in the context of this 
work, the term loss of identity does not refer to a parton 
change of color or phase, which are accounted for already 
in the description of the underlying QCD scattering pro- 
cess, neither to a change of flavor which would be me- 
diated by the weak interaction and thus irrelevant for 
the time scales involved during the QCD plasma phase. 
This term is rather related to the onset of a hadroniza- 
tion mechanism that happens during the interaction of 
a fast (hard) parton with soft ones from the medium in 
such a way as to produce a hadron by means of recombi- 
nation. This process has been referred to as shower and 
thermal parton recombination, in Ref. (22J. In this way, 
the loss of identity is related to the fact that when this 
kind of partons form a hadron, the energy loss should 
stop being described in terms of parton degrees of free- 
dom and start being described in terms of hadron degrees 
of freedom, thereby effectively producing the parton to 
disappear from the description. 

Typical time scales involved in hard-soft parton recom- 
bination are of order T recomb ~ 1.5 fm, despite the low mo- 
mentum transfers involved, since, as argued for instance 
in Ref. [23|, this recombination needs not be local and 
it can be mediated by a QCD string. Such time scale is 
well within the life-time of the QCD medium, for central 
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collisions and the largest nuclei, where it is estimated to 
be of order 5 fm. 

When the medium length shortens and falls below 
Trccomb, hadronization is more likely to happen outside 
the medium. This means in particular that for periph- 
eral collisions or collisions of smaller systems, the energy 
loss description in terms of partonic degrees of freedom 
is appropriate. This is accounted for in our description 
when we take r^ 1 > L, for which the energy loss with 
and without the use of the parameter rj coincide. Indeed, 
a change in the energy loss of intermediate momentum 
hadrons should exist as a function of centrality and sys- 
tem size. In order to quantify such change it is necessary 
to quantitatively estimate the energy loss of a hadron 
within a QCD medium. This calculation is for the mo- 
ment outside the scope of the present work. 

In the context of recombination, the use of the parame- 
ter rj should be introduced into a statistical scenario that 
also incorporates the evolution of the colliding system 
with energy density, as well as using a realistic geom- 
etry, including adequate probability profiles to produce 
jets [24| and the effects of an expanding medium. All this 



is for future. 
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V. APPENDIX 

The functions defined in Eq. (|21|) are explicitly given 

by 
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, dQ 
1 



sin[(w v • q + irj)L/2v] 



8qvr] 



u> — v • q + ir] 



icosh(7yL/w) sgn(Z)sgn(Z')Cz[(cj — Iqv — l'ir])L/v} 

i,i'=±i 



+ 2i 



arctan 



qv + co 



+ arctan 



qv — uj 
V 



+ smh(r]L/v) sgn(l)Si[(oj — Iqv — l'ir])L/v] 



,U'=±i 



Ji 



47T 



1 

8qv 



sinhw — v • q + in)L/2v] 

— (u - v • q) 

uj — v • q + irj 

cosh(r/L/v) sgn(l)Ci[(oj — Iqv — l'irj)L/v] 



K i,i'=±i 



+ In 



(uj + qv) 2 + rj 2 



(uj — qv) 2 + rj 2 _ 
smh(-qL/v)\ ^ sgn(l)sgn(l')Si[(uj - Iqv - l'ir\)Ljv\ 



J2 



f dn 

J 4^ 

1 

8qv 



sin[(w v • q + i-q)L/2v] 



(uj-v q) 2 



uj — v • q + ir] 

rjcosh.(rjL/v) \ i ^ sgn(Z)sgn(Z')Ci[(w — Iqv — l'ir])L/v] 



arctan 



qv + ui 
V 



=±i 



+ arctan 



qv — ui 
V 



4qv 
V 



+ rjsmh(r)L/v) I sgn(l)Si[(uj — Iqv — l'ir])L/v} 

\m'=±i 

— I — J cos(Luj/v) sin(Lq) 



(27) 



I 

where sng is the sign function and Ci and are the pearance, the above functions are all real for real values 
cosine and sine integrals, respectively. Despite their ap- of u>, q, v, L and rj. 
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